EFFECTIVE FREQUENCY TECHNIQUE 


C. Laurence Korb 

Laboratory for Atmospheres, Code 912 
NASA/Goddard Space Flight Center 
Greenbelt, MD 2077 1 

and 


Chi Y. Weng 

Science Systems and Applications, Inc. 
Seabrook, MD 20706 



Abstract 


An effective monochromatic frequency technique is described to represent the effects of 
finite spectral bandwidth for active and passive measurements centered on an absorption 
line, a trough region, or a slowly varying spectral feature. For Gaussian and rectangular 
laser line shapes, the effective frequency is shown to have a simple form which depends 
only on the instrumental line shape and bandwidth and not on the absorption line profile. 
The technique yields accuracies better than 0. 1 % for bandwidths less than 0.2 times the 
atmospheric line width. 



I. Introduction 

The effects of finite bandwidth distort a spectral measurement by broadening the line and 
decreasing the apparent absorption at line center. This can produce large errors in the 
measured spectrum . 1 For laboratory measurements, deconvolution techniques may be 
used to recover the true spectrum . 2 The analysis of range resolved lidar measurements of 
differential absorption in the atmosphere is more complex since the spectral line shape 
changes with altitude and the intervening atmosphere affects the measured differential 
absorption at the end point of the measurement . 3,4 A method of solving for the 
differential absorption has been briefly discussed using a complex iteration procedure . 3,5 

In this paper, we describe a fundamentally new approach to analyze the problem of 
radiative transfer and atmospheric transmission with finite spectral bandwidth. We show that 
an equivalent monochromatic frequency may be used to represent the effects of finite 
bandwidth. In the vicinity of line center, the equivalent frequency can be expressed by a 
simple analytical form that depends on the characteristics of the instrumental distortion 
function rather than the nature of the spectral line profile and the particular characteristics of 
the path being used. The technique is applicable to radiative transfer calculations as well as to 
a large class of active and passive spectroscopic experiments, in both the laboratory and the 
atmosphere. An equivalent frequency not only greatly simplifies the computation of finite 
bandwidth effects for forward problems in radiative transfer but also allows the inverse 
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problem of the remote measurement of species concentration or atmospheric properties to be 
solved using simplified monochromatic methods. 

An integrated form of the mean value theorem provides the theoretical basis for this 
technique 6 . It can be shown for a finite bandwidth measurement that an equivalent frequency 
exists at some point within the instrumental bandwidth if the instrumental line profile is a 
continuous function and does not change sign. The location of the equivalent frequency is not 
defined, however. In this paper, we present general methods for calculating the equivalent 
frequency and show that the results are applicable to the whole range of spectroscopic and 
atmospheric line profiles. 

II. Theory 

For the case of finite bandwidth, the two way transmission of the atmosphere is related to the 
absorption coefficient K(d,x), the absorption per unit path as 3 

T u (0,/?) = J /i(t;)exp|-2| o ^(t),x)<ijcJii) (1) 

where the integral of the instrumental line shape h(u) is normalized to have unit value. In the 
vicinity of an absorption line centered at v 0 , the absorption coefficient may be expanded as 

K('o,x)=a 1 (x)-a 2 (x) (o-d 0 ) 2 + o ((v-v 0 ) 4 ) (2) 

where we explicitly allow for a dependence of the absorption coefficient on distance as, for 
example, would be encountered in the atmosphere. The quadratic expression of Eq. (2) is valid 
in general since a Taylor series expansion of any function about a maximum or minimum 
yields an equation of this functional form whereas dK/d'ol "0 o =0 and the linear term in the 
expansion is zero. Eq. (2) provides a good approximation to the absorption coefficient for 
bandwidths which are small compared to the absorption line width. For a collision-broadened 
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or Doppler profile, Eq. (2) gives the absorption coefficient to accuracies better than 1% 
for distances up to 0.3 half-widths from the line center. Since Eq. (2) is general, it can be 
used to represent the absorption coefficient near line center for the Voigt line profile or 
even for a complex line shape such as would result by integrating Voigt profiles over an 
atmospheric path with varying pressure and temperature. Eq. (2) can also be used to 
represent the absorption in trough regions, the region of minimum absorption between 
two absorption lines. 7 

To evaluate the transmission for a measurement centered on an absorption line, an 
explicit form for the observational line shape is required. We consider first the case of a 
Gaussian profile, e.g.. 


h(t*= 


_ 2 _ 

Av 



exp (- 1 n 2 [('D-u 0 )/AW2] 2 ) 


(3) 


It follows from Eqs. (1), (2), and (3) that the transmission is given as 

t u ( 0,R)=(^) exp(-2j o a 1 (x)dx)j_exp{-('u-'0 0 ) 2 [b-2}2i 2 (x)dx]}d'u (4) 


where b=4 In2/At> 2 . Eq. (4) can be integrated and for 

2 fa 2 (x)dx/b«l, (5) 

Jo 

the result can be expressed as 

t u (0,R)=exp{-2 f* [ a,(x)-a 2 (x)/(2b)(l+ [ a 2 (x)dx/b)]dx} (6) 

JO JO 
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Eq. (5) is equivalent to the condition that the change in the optical depth between the 
center and half-width point of the measurement bandpass is small with respect to ln2. It 
then follows from Eqs. (1) and (2) that Eq. (6) is equivalent to 

t v (0,R)=exp(-2j ( A' ( u, x)dx), 


where 


t) e = u ±Ad/(2V21«2 )( 


1 + 


At ) 2 
81n2 



(7) 


Thus, the transmission for a measurement centered on an atmospheric line or in a trough 
region can be represented by an equivalent monochromatic measurement at the frequency 
\) e for bandwidths which are small compared to the atmospheric line width. We note 
that the first term for the equivalent frequency depends only on the characteristics of the 
instrumental distortion function. 

An effective frequency can also be used to analyze finite bandwidth measurements for a 
measurement with a rectangular line shape. For this case, the line shape is 
h(u)= 1/Av | \)-\) 0 1 <_Av/2 

= 0 | i)-o 0 1 >At)/2 (8) 

— fAu/2 . 

and t v (0,R)=(1/Ai)) exp[-2 (a,(x)-a 2 (x)(n-'O 0 ) 2 )dx]d\) 

J-Av/2 JO 

We note that a rectangular line shape provides a first order approximation to the average 
spectral output of a multimode laser with a frequency jitter of the order of the mode 
separation. From Eqs. (1), (2), (8) and an analysis similar to that given for Eqs. (4)-(7), it 
can be shown that for a rectangular line shape 

u=u 0 ±Au/(2V3) (9) 
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Eqs. (7) and (9) show that for a Gaussian or rectangluar instrumental line shape, the 
transmission over the atmospheric path from range 0 to R can be represented simply by 
the transmission at the effective frequency u e . We note that Eqs. (7) and (9) are valid 
for any atmospheric path. 

We can also utilize the equivalent frequency for analysis of finite bandwidth effects at a 
reference frequency. Since a reference frequency is generally chosen either in the far 
wing of a line or on the edge of an absorption band where the absorption changes slowly 
with frequency, the absorption coefficient may be represented as a linear function of 
frequency for a moderately narrow bandwidth. That is, 

K(u,x)=a(x) + b(x)(D-D r ). (10) 

For the case where the variation in the atmospheric optical depth is small over the line 
shape, i.e., 

2('U-\) r ) Jl>(x) dx«l, (11) 

it follows that the corresponding term in the transmission in Eq. (1) can be expanded 
which yields 

t r (0.R)=exp(-2 fa(x)dx) f h 1) [l-2(t>-'U r ) f b(x)dx]du (12) 

JO J-oo Jo 

For a line shape which is symmetric about i) r , h (u) Cu-\) r ) is an odd function (e.g., with 
zero integral) and Eq. (12) reduces to 

t r (0,R)=exp(-2 \K (u r ,x)dx) (13) 

and o e =D r 

Thus, the transmission in the region of the reference frequency can be represented by an 
equivalent monochromatic measurement at the center of the bandpass for any 
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instrumental line shape which is symmetric providing the absorption coefficient is a 
linear function of frequency. 

III. Application to Differential Absorption 

Consider a two-frequency lidar measurement of the resonant absorption for an 
atmospheric species (e.g., C0 2 or H 2 0) or property (e.g., temperature or pressure) where 
the measurement at each frequency is made with a finite bandwidth laser. The 
measurement at the on-line frequency v is typically centered on an absorption line of the 
desired atmospheric species. A reference frequency u r is then selected at a nearby 
location where the resonant absorption is weak compared to that at the measurement 
frequency. The two frequencies are chosen sufficiently close, however, so that the 
atmospheric backscatter and the attenuation due to both scattering and continuum 
absorption processes is essentially identical at the two frequencies. The ratio of the 
signal returns backscattered from the atmosphere at two adjacent frequencies at ranges R, 
and R 2 essentially eliminates the scattering and absorption properties of the atmosphere, 
except for those of the resonant absorption effect to be measured. For elastic scattering 
and homogeneous scattering and absorption over a range gate of width AR, it follows 
from the lidar equation 8 that the ratio R of the energy at the on-line and reference 
frequencies from two successive range gates is given as 

t.(0,R,)/z(0,R,) 

where 

R £■,(/?,)/&(&) 
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and where x u (0,R) and t r (0,R) are the transmission along the two-way path path 
from 0 to R for laser bandwidths centered at the on-line frequencies o and the reference 
frequency u r , respectively, and e (R) is the energy received from an element of the 
atmosphere at range R with thickness AR, and R 2 =R,+AR. 

* 

The inclusion of finite laser bandwidth effects and non-zero resonant absorption 
at the reference frequency complicates the analysis of Eq. (14) considerably since the 
differential absorption coefficient, K(\>) - K(u r ) is not a directly measured or easily 
derived quantity. In addition, the differential ranging measurement given by Eq. (14) 
not only has absorption information in the desired range element of width AR centered at 
R=(R,+ R 2 )/2, but also contains absorption information from the region 0 to R,. 


If an effective frequency is used to represent the absorption at the reference frequency 
(i.e., Eq. (13)), it follows that a differential absorption measurement with differential 
ranging and finite bandwidth at both the on-line reference frequencies can be represented 


f“ A(u)ex pf-2 f*fr(t>,*) - K(v„x)]dx\dv 

* = “ Hk f- 09 

J h(v)exp<-2jJ.K(v,x)- K(Vr,x)]dx ■ dv 


where h(o) is the laser line shape for the on-line measurement. Eq. (15) has the form of a 
differential absorption experiment. We note that Eq. (15) is valid even for the case of 
strong absorption in the region of u r since no assumptions on the overall level of 
absorption were made. 
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An effective frequency can also be used to greatly simplify the analysis of finite 
bandwidth measurement at the on-line frequency. It follows from Eq. (15) that a 
differentia] absorption measurement with differential ranging and finite bandwidth can be 
expressed in terms of the effective frequencies v e and D r as 

R = exp {-2 f* 2 [ K(u,x)-K( D r ,x)]dx). (16) 

Thus, the effective frequency method reduces the problem of measurements with finite 
bandwidth at the on-line and reference measurements from a relatively untractable form, 
where the measured quantity is a complex function and depends on the absorption over 
the entire atmosphere path, to a simple measurement of the differential absorption 
coefficient, K (t> e ,x) -K(\) r ,x), over the desired measurement interval from R,to R 2 . 

IV. Comparison with Simulation 

To evaluate the application of the effective frequency technique, we consider the case of 
generalized absorption in the atmosphere. The absorption coefficient of a gas species in 
the atmosphere at frequency t) is 3 

K(D)=qnS(T)f(D-U 0 ), (17) 

where q is the molecular mixing ratio of the species, n is the number density at pressure p 
and temperature T, S (T) is the line strength, and f (v-v 0 ) is the line shape at frequency u 
for a line centered at 0) o . For atmospheric measurements where the effects of both 
collision and Doppler broadening contribute significantly to the absorption, the line shape 
is given by the Voigt profile 9 which is a convolution of independent Lorentz and Doppler 
profiles. It is given as 
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(18) 


where 


f(U-D 0 )= 


fa 

n 


f. 


exp(-y 2 )Jy 

a 2 +(C-y) 2 ’ 


f=( ln2/7r) 1/2 b d , 
a=(ln2) 1/2 b c /b d , 

Hln2) 1/2 (D-v 0 )/bd (19) 

and b c and bd are the collision broadened and Doppler broadened line half-widths at half- 
height, respectively. These are defined as 

bc(p,T)= b° c (To)(p/po)(To/T) m , 

b d (T)=[2kTl n2/(M 0 c 2 )] 1/2 u 0 , (20) 

where Mo is the mass of the absorbing species, c is the velocity of light, po is standard 
pressure, b" (T) is the collision-broadened line half-width at a pressure of one 
atmosphere and the exponent of the temperature dependence of the collision-broadened 
half- width m is predicted to be 1/2 from kinetic theory. 


To evaluate the error in the determination of the differential absorption coefficient using 
an effective frequency, Eq. (15) was evaluated numerically for rectangular and Gaussian 
laser line shapes of various widths using Eqs. (17) and (18). The atmospheric model used 
for evaluating Eq. (17) is a multilayer, one-dimensional model with up to 50 
homogeneous layers. The calculations simulated measurements for a two-way 
atmospheric path from infinity to various altitudes for a total optical depth of two at line 
center. The altitude of the end point of the measurement is represented by the pressure at 
that altitude as given by the parameter “a” which is proportional to pressure. The 
calculations thus begin with a Doppler profile at high altitude, a=0, and include all values 
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of “a” for the Voigt line profile up to and including the value of ‘a” at the end point of 
measurement. 

Figures 1 A and IB give the error between this detailed evaluation of the differential 
absorption coefficient and the analysis using a monochromatic calculation at the effective 
frequency for the cases of Gaussian and rectangular laser line shapes, respectively. The 
errors are shown as a function of the ratio of the laser bandwidth, the full width at half 
height, to the half- width at half-height of the Voigt line profile at the end point of the 
measurement, b v and a. For values of Au/b u less than 0.5, the error in using an effective 
frequency is less than 0.5% for a rectangular line shape and less than 1% for a Gaussian. 
Accuracy at the 0.1% level can be obtained for values of Au/b u less than 0.2 for a 
rectangular line shape or at the 0.2% level for a Gaussian.. 

V. Conclusion 

In conclusion, we have shown that an effective monochromatic frequency can be used to 
represent the effects of finite instrumental bandwidth for spectral measurements centered 
on an absorption line, a trough region, or a slowly varying spectral feature. The 
technique is applicable to both integrated path and differential ranging experiments and 
the resulting frequency is independent of the atmospheric line shape and path, to first 
order. The application of an effective frequency technique greatly simplifies the analysis 
for a large number of spectroscopic problems. Analyses similar to those we have 
presented for Gaussian and rectangular line shapes can be used to derive the effective 
frequency for other instrumental line shapes. Complex line shapes, including multimode 
structure, can also be represented by a set of effective frequencies which correspond to a 
superposition of elementary forms. 
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Figure Captions 


Figure 1A and IB. The percentage error in the differential absorption coefficient, as a 
function of the laser bandwidth to the half- width at half-height of the Voigt line profile, 
using a monochromatic calculation at an effective frequency to represent the effects of 
finite laser bandwidth. Figures 1 A and IB are for the cases of Gaussian and rectangular 
laser line shapes, respectively, for a downward viewing experiment from space to an 
altitude (pressure) in the atmosphere defined by the Voigt parameter a=(b° c p/bd)Vln2. 
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POPULAR SUMMARY 


This paper describes a fundamentally new and simple approach to calculating the percent 
of energy transmitted over any atmospheric path or through a gas cell. These calculations 
are important for determining the amount of pollution in the atmosphere or for 
determining the amount of atmospheric heating due to ‘the greenhouse effect’ . Normally, 
calculations such as these are complex and require a computer to do thousands of 
computations to represent the effects of each color of the spectrum that is involved in the 
process. These calculations usually also depend on many other characteristics such as the 
particular atmospheric path, the characteristics of the absorbing gas, and the 
characteristics of the instrument used to do the measurement. 

We have shown that rather than having to do these calculations for each color of the 
spectrum that is involved, we only have to do the calculations at a single effective color 
that represents the average effect of all the colors involved in the process. We have 
shown that this average effective color has a simple form that depends only on the nature 
of the averaging or measuring process and not on the characteristics of the atmospheric 
path or the absorbing gas. 

This technique is applicable to a large class of calculations and measurements both in the 
laboratory and the atmosphere. This method allows accuracies in the calculation or 
measurement of the percent of energy transmitted by a gas to be as high as 99.9%. 
Simplifying the calculation and analysis of measurements reduces the need for thousands 
of computations to a single computation that a high school student could do with a good 
calculator. 


